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Abstract. We investigate the variability of the QPO frequency in 4U 1608- 
52 on very short time scales of a few seconds. We detect changes in frequency 
as fast as 0.5Hz/s, which is more than ten times the fastest change reported 
for this object. Using a structure function analysis, we are able to detect 
\ variability on time scales as short as ~ 4 s. We do not detect any time 

. asymmetry in the QPO-frequency history. These results have important 

consequences on the measurement of the QPO coherence. 
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Fourier analysis has shown that the X-ray emission of low-mass X-ray binaries 
(LMXBs) often presents characteristic quasi-periodic-oscillation (QPO) peaks at 
^ ' frequencies up to 1 kHz (van der Klis 1989, 2000). The kHz QPOs are of par- 
• ticular interest, as they are a probable signature of physical mechanisms taking 
^ place in the inner parts of the accretion disk, where strong-field effects of gravity 
become important. The physical origin of the QPOs is still a matter of debate, 
and a detailed understanding of the phenomenon is therefore necessary. Berger 
^ , et al. (1996) analyzed the variation of the QPO frequency in a Rossi- XTE ob- 
■ ■ ' servation and observed a 40 Hz change in frequency in about 1000 s. They also 
found evidence for a high QPO coherence of the order of 200. Building upon these 
results. Barret et al. (2005) (hereafter Paper I) have recently made important 
progress in the determination of the width of the kHz QPO in 4U 1608-52. In 
Paper I they found evidence that the variability of the QPO frequency could be 
quite significant, even on short time scales. In this work, we analyze the temporal 
properties of the QPO frequency in 4U 1608-52, focussing on the shortest possible 
time scales allowed by the signal. 
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2 Data 



We use observations of 4U 1608-52 from the Rossi-XTE public archive. These 

data have aheady been analyzed in Berger et al. (1996) and Mendez et al. (1998a, 
1998b, 1999), as well as in Paper I. We selected ~ 1000 seconds from an observation 
segment of 4U 1608-52 performed on March 3rd, 1996 20:54. This data set has 
been chosen because the QPO is particularly strong and highly coherent, and its 
frequency has been found in Paper I to remain in a narrow frequency range. 

3 Analysis methods 

3.1 Rayleigh's test as an estimator of the PSD 

We consider the individual photons, without building a light curve. The observed 
signal is therefore a sum of Dirac 6 functions: 



where the ti's are the time tags of the individual photons. The Fourier transform 
and the powcr-spcctnmi density (PSD) of this signal can be easily obtained: 



FT(5) = ^exp(27rz/(.fi) ^ PSD(S') = || ^exp(27rz/(.ii)||2 (3.2) 



This is exactly the expression of Rayleigh's test (e.g.. Fisher 1993). We shall use 
this expression to estimate the PSD of the light curve of 4U 1608-52, in order to 
avoid any loss of information due to the binning of the photons. 

3.2 QPO-frequency determination 

We determine the variation of the QPO frequency using a simple method: We 
calculate the PSD over time intervals as short as a few seconds, and smooth 
it using a Gaussian; we then search the maximum of the smoothed PSD. This 
method has been validated using simulations of a QPO at a known frequency with 
similar amplitude and width (sec Paper I) . Wc find that the accuracy of the QPO 
frequency determination depends little on the width of the Gaussian, provided 
that it is significantly wider than the natural resolution of the Fourier transform, 
i.e. Ajy ^ l/T, where T is the duration of the time interval. We used a Gaussian 
with a FWHM of 3 Hz for time intervals from 1 to 10 s. 

Statistical fluctuations may sometimes produce spurious peaks which exceed 
the strength of the QPO. To eliminate these peaks, wc first calculate the PSDs 
from 1/ = 100 up to = 2000 Hz over 10s time bins. With such time intervals the 
measured QPO frequency always fall within a few Hz from the curve determined in 
Paper I. When we shorten the time bin, spurious peaks appear. To get rid of these 
peaks, we first restrict the range of investigated frequencies to i/ = 810 — 850 Hz, 
and reject QPO-frequency measurements that fall farther than 5 Hz away from the 
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Fig. 1. (a) Evolution of the QPO frequency calculated every 1 s over 10s bins, (b) Zoom 
on the QPO-frequency evolution in 10s (black curve) and 4s (grey curve) bins. In each 
panel the average uncertainties are shown in the upper right corner. 

lOs-bin frequency. Using this procedure we eliminate about 2% of the peaks when 
using 4s bins. This figure increases to 11% for 2s bins, and to 25% for Is bins. 

4 Results 

4.1 QPO frequency as a function of time 

Fig. ^ shows the evolution of the QPO frequency as a function of time, i'Qpo(i), 
for the full segment of data analyzed here. The QPO frequency has been de- 
termined every 1 second using 10 s time bins. Very rapid changes of frequency 
as large as 0.5Hz/s are observed, both towards higher and lower frequencies. A 
Kolmogorov-Smirnov test on the distributions of the positive and negative fre- 
quency changes does not show any evidence of time asymmetry. Fig. zooms 
in on a portion of the data, on which we superimpose VQPoit) calculated in 10s 
and 4s bins. While the curves do follow each other quite well, i^Qpo{t) with 4s 
bins show much more prominent structures that it seems difficult to ascribe to the 
increased uncertainties only. 

4.2 Time-series properties 

To investigate and quantify the variability of the QPO frequency, we compute the 
structure functions (SF) of t'QpoCi) for the data and for simulations. The SF of 
a time series (e.g., Paltani 1999) measures the amount of variability as a function 
of the time scale r. It is flat for a white-noise-dominated time series. The SF 
increases with t as long as the time series does contain some variability at the 
corresponding time scale r. Fig. shows the SFs of i'Qpo(i) for different time 
bins. For comparison Fig. shows the same SFs for the simulations. These SFs 
are dominated by the uncertainty in the frequency determination, which shows 
that our method does not produce any artefact. For the 4U 1608-52 data, the 
horizontal part is present in the SF of vqpo {t) calculated in 1 s bins only. Changes 
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Fig. 2. SFs of fQpo(t) for the observed (a) and simulated (b) data for different bin 
sizes. The bin size increases from top to bottom. 

in the QPO frequency on time scales as short as 4 s are therefore detected. The 
differences in the SFs are the result of the smoothing and of the decrease in the 
uncertainty introduced by the use of larger bins. 

5 Conclusion 

We observe the presence of an erratic variability of the QPO frequency in 4U 1608- 

52, with evidence of changes in QPO frequency as large as 0.5 Hz/s. This is already 
more than ten times faster than previously reported by Berger et al. (1996). 
Faster frequency changes seem present in z>'Qpo(t) calculated with shorter time 
bins. The signal is strong enough, and the method developed here sufficiently 
accurate, that a minimum variability time scale of ~ 4s is detected. We did not 
find any evidence for asymmetry in UQPo{t). These rapid and erratic frequency 
changes have consequences for the measurement of QPO coherence. A blurring 
of the QPO profile will indeed take place if too long time bins are used, leading 
to an underestimation of the true QPO coherence. As discussed in Paper I, the 
high coherence of the signal places very strong constraints on QPO models. It is 
therefore important to measure it on very short time scales. 
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